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Predicting the geographic distribution of Lutzomyia longipalpis 
(Diptera: Psychodidae) and visceral leishmaniasis in the state 
of Mato Grosso do Sul, Brazil 
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To understand the geographic distribution of visceral leishmaniasis (VL) in the state of Mato Grosso do Sul (MS), 
Brazil, both the climatic niches o/Lutzomyia longipalpis and VL cases were analysed. Distributional data were ob- 
tained from 55 of the 79 counties of MS between 2003-2012. Ecological niche models (ENM) o/Lu. longipalpis and 
VL cases were produced using the maximum entropy algorithm based on eight climatic variables. Lu. longipalpis 
showed a wide distribution in MS. The highest climatic suitability for Lu. longipalpis was observed in southern MS. 
Temperature seasonality and annual mean precipitation were the variables that most influenced these models. Two 
areas of high climatic suitability for the occurrence ofVL cases were predicted: one near Aquidauana and another 
encompassing several municipalities in the southeast region of MS. As expected, a large overlap between the models 
for Lu. longipalpis and VL cases was detected. Northern and northM'estern areas of MS were suitable for the oc- 
currence of cases, but did not show high climatic suitability for Lu. longipalpis. ENM of vectors and human cases 
provided a greater understanding of the geographic distribution ofVL in MS, which can be applied to the develop- 
ment of future surveillance strategies. 
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Visceral leishmaniasis (VL) is a severe systemic dis- 
ease resulting from infection with Leishmania infantum 
(syn. Leishmania chagasi). In the Americas, domestic 
dogs are the main reservoirs of the parasite, which is 
transmitted mainly by Lutzomyia longipalpis (Diptera: 
Psychodidae) (Gontijo & Melo 2004), a synanthropic 
sandfly species with a wide geographic distribution in 
Brazil (Souza et al. 2009). Between 2003-2007, approxi- 
mately 3,500 VL cases were reported per year in Bra- 
zil, which corresponds to 95% of all new cases of VL in 
the Americas (Alvar et al. 2012). Control measures di- 
rected against sandflies have been unsuccessful and the 
early detection and treatment of human cases remains 
the most important strategy for reducing case fatalities 
(Romero & Boelaert 2010). 

In the state of Mato Grosso do Sul (MS), VL was 
restricted to the municipalities of Corumba and Ladario 
until 1993 (Antonialli et al. 2007), but its presence is 
currently registered in 41 municipalities. An increase 
in the number of cases diagnosed in urban centres was 
observed, reflecting both disease expansion and ur- 
banisation in the state as well as representing a threat 
to public health, considering the mean disease lethality 
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of 7.8% (Oliveira et al. 2006b). The Notifiable Diseases 
Information System recorded 1,919 cases of VL in MS 
between 2003-2011, which represents 15% of all cases 
recorded in Central-West Brazil. Most cases occurred in 
the municipalities of Aquidauana, Campo Grande and 
Tres Lagoas, which are currently classified as areas of 
intense transmission of VL, where Lu. longipalpis has 
also been detected (Silva et al. 2007, MS/SVS 2013). 

The geographic distribution of Lu. longipalpis ex- 
tends from Mexico to Argentina (Deane & Deane 1962, 
Sherlock 1996, Souza et al. 2009, Salomon et al. 2011). 
This species occurs in areas with different climatic con- 
ditions, ranging from the semi-arid areas in the Caatin- 
ga biome to wetter areas in other Brazilian biomes and 
presents great adaptability to different habitats and en- 
vironmental changes (Ximenes et al. 2000, Aguiar et al. 
2003, Missawa & Lima 2007). In MS, Lu. longipalpis is 
the primary species captured in anthropic environments 
(Oliveira et al. 2006b, Almeida et al. 2010a, b) and is the 
principal vector of Z. infantum. Transmission, however, 
can occur by Lutzomyia cruzi in Corumba in northwest 
MS (Santos et al. 1998). 

Ecological niche models (ENM) have been used to 
predict the occurrence areas of sandflies and leishma- 
niasis in Brazil. Peterson and Shaw (2003) predicted 
climate change effects on Nyssomyia whitmani distri- 
bution; specifically, under climate change scenarios, an 
extension of the geographic range oiNy. whitmani was 
predicted in Southeast Brazil, where cutaneous leishma- 
niasis appears to be re-emerging. Peterson et al. (2004) 
demonstrated that elevation and vegetation significantly 
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influence the distribution of sandflies in the state of Sao 
Paulo (SP). In the state of Bahia, predictive models es- 
timated a higher incidence of VL in hot and dry areas 
(Caatinga), whereas areas of Cerrado and the Atlantic 
Forest were considered to be less suitable for disease 
transmission (Nieto et al. 2006). Through logistic re- 
gression, Zeilhofer et al. (2008) demonstrated that the 
spatial distribution of Ny. whitmani is associated with 
deforestation and with lower rates of gross domestic 
product in the state of Mato Grosso (MT). 

Given the proven usefulness of spatial modelling 
techniques in understanding the distribution of sandflies 
and leishmaniasis and the fact that the predictive ana- 
lysis of the geographic distribution of Lu. longipalpis 
and human cases of VL can aid in the development of 
monitoring and disease control programs in MS, the 
present study aims to analyse the climatic conditions of 
the occurrence areas of Lu. longipalpis and VL cases 
using ecological niche modelling. 

MATERIALS AND METHODS 

MS is located in the Central- West Region of Bra- 
zil, with an area of 358,200 km^. It borders the states 
of Minas Gerais, SP and Parana to the east and MT and 
Goias to the north. Moreover, MS has international bor- 
ders with Bolivia and Paraguay. The state has 78 coun- 
ties, divided into 11 micro-regions, with a population of 
approximately two million inhabitants (IBGE/MPOG 
2010). The Sedimentary Plateau Basin of the Parana 
River and the South Pantanal Plains are the major geo- 
logical units in MS. The biomes present in the state are 
the Cerrado, Pantanal and Atlantic Forests. 

The climate of MS, according to the Koppen classi- 
fication, varies between the Aw (tropical humid climate 
with a rainy summer and a dry winter) and the Cfa (hu- 
mid mesothermal climate with no dry season) subtypes, 
with temperatures during the warmest month reaching 
over 22°C and with low precipitation during its driest 
month. In MS, there is predominance of undulated el- 
evation and of Oxisols (IBGE 1990). 

Distributional data for Lu. longipalpis were obtained 
from entomological surveys carried out by the Sector of 
Entomology of the Vector Control Coordination of the 
State Department of Health and from entomology labo- 
ratories of the respective municipalities. The research 
was conducted between 2003-2012 in urbanised areas of 
55 municipalities of MS, corresponding to 70% of the 
counties in this state. Samples were collected with CDC 
light traps that were installed at dusk inside houses and 
in domestic animal shelters (peridomiciles) and collec- 
tions took place over a period of 12 h for three consecu- 
tive nights. In 23 of the 55 municipalities (42%), a survey 
of sandflies was performed via entomological monitor- 
ing over two years following the methods described in 
Almeida et al. (2010b). The number of points sampled 
varied among municipalities depending on the urban 
total area (between 4-12 sampling points), with capture 
effort ranging from 96-1,728 traps by municipality. For 
other municipalities, entomological surveys were con- 
ducted with an average capture effort of 50 traps. For 
identification of sandfly species, specific keys were used 



(Galati 2003). Sandfly specimens were deposited in the 
Regional Laboratory of Entomology of Dourados, MS. 
The presence oiLu. longipalpis was recorded in 118 sites 
and their georeferencing had a precision of 0.0001°. 

The present study included records of VL cases re- 
ported by local health services in the municipalities of 
MS during 2003-2012. The cases were georeferenced 
according to the coordinates of the municipality admin- 
istrative centre with a reasonable degree of confidence 
(i.e., with an uncertainty of < 5 km, to a precision of ± 
0.01°). All records were georeferenced based on a consul- 
tation of Brazilian Institute of Geography and Statistic 
(IBGE) (ibge.gov.br/).To characterise climatic variation 
across MS, eight climatic variables were used: annual 
mean temperature, diurnal temperature range, tempera- 
ture seasonality, maximum temperature in the warmest 
month, minimum temperature in the coldest month, an- 
nual precipitation, precipitation in the wettest month and 
precipitation in the driest month. These variables were 
obtained from the WorldClim project (worldclim.org) 
and were developed via the interpolation of mean month- 
ly climatic data from meteorological stations over 30-50 
years (1950-2000), depending on data availability at sta- 
tions (Hijmans et al. 2005). These variables were selected 
to avoid the confounding effects of calibrating models in 
an overly dimensional environmental space (Peterson & 
Nakazawa 2008). The climatic database used in the an- 
alysis presented a spatial resolution of 5 x 5 km per pixel. 

Models were produced using Maxent v.3.2.1 (Phil- 
lips et al. 2006), applying default parameters with the 
exception that we chose a random seed, with 10 replicate 
analyses based on bootstrap subsampling. We used the 
median output grids as the best hypothesis of potential 
range and imported them into ArcView v.3.3 for assess- 
ment and analysis. Distributional data for Lu. longipalpis 
and human cases of VL were separated into two sets: one 
for model calibration (75% of points) and the other for 
model evaluation {25% of points). For Lu. longipalpis, 
several sites were near one another and the model was 
produced by using a reduced number of spatially unique 
points, considering the spatial resolution of our data. 

The quality of the models generated was evaluated 
using the receiver operating characteristics (ROC) curve, 
which correlates the sensitivity and specificity of mod- 
el predictions across multiple thresholds (Phillips et al. 
2006). The area under the ROC curve (AUC) provides 
a measure of model performance across many thresh- 
olds. AUC values close to one indicate high performance, 
whereas readings below 0.5 indicate the poor perfor- 
mance of a model (Elith et al. 2006). Model accuracy was 
assessed by examining the omission rates associated with 
test points (Anderson et al. 2002) and Maxent's jackknife 
test was used to identify variables that most influence 
model predictions. The contribution of each variable to 
the model was estimated through the permutation impor- 
tance provided by Maxent (Phillips et al. 2006). To test 
model significance, we compared the predictive success 
of the models against null expectations using a cumu- 
lative binomial test. In particular, we assessed whether 
each test point fell in the areas that were identified by the 
model as suitable and compared this success rate with 
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overall proportions of pixels identified as suitable or un- 
suitable for that species. 

Finally, to analyse the overlapping area between the 
models of Lii. longipalpis and VL cases, the raw Maxent 
outputs were converted into binary maps showing the 
presence or absence based on a threshold of 0.60 (values 
above this threshold were considered presence and they 
represented higher climate suitability). Then, the over- 
lapping area was calculated using the ArcView Spatial 
Analyst extension. 

RESULTS 

Lu. longipalpis occurred in 37 (67%) of the 55 mu- 
nicipalities surveyed. The predictive model indicated 
wide potential distribution oiLu. longipalpis in MS. The 
highest climatic suitability for Lu. longipalpis was ob- 
served in the south of the state (Fig. 1). Lu. longipalpis 
was not registered in areas with high climatic suitability 
for vector occurrence in several municipalities of MS, 
such as Navirai, Itaquira and Juti. Temperature season- 
ality and annual mean precipitation were the variables 
that most influenced the distribution of Lu. longipalpis 
(Maxent's permutation importance = 79.7 and 4.9, re- 
spectively). The model revealed 0% error of omission 
and AUC = 0.992 and was statistically significant (p < 
0.01) as indicated by the binomial test. 

Two areas of high climatic suitability for the occur- 
rence of VL cases were predicted: one near Aquidauana, 
Dois Irmaos do Buriti and Anastacio and another en- 
compassing several municipalities in the southeast re- 
gion of the state, such as Nova Alvorada do Sul, An- 
gelica and Jatei (Fig. 1). Temperature seasonality and 
precipitation of the driest month were the variables that 
most influenced models of VL cases (Maxent's permuta- 
tion importance = 67.4 and 7.3, respectively). The model 
showed 0% error of omission and AUC = 0.991 and was 
also statistically significant (p < 0.01). 




Fig 1: ecological niche models projected for the potential distribu- 
tion of Lutzomyia longipalpis (left) and human cases of visceral leish- 
maniasis (right). Areas identified as suitable based on the records of 
the occurrence of Lu. longipalpis are shown on a red scale: light red 
(low suitability) to dark red (high suitability), while areas identified 
as suitable based only on VL cases are presented on a green scale: 
light green (low suitability) to dark green (high suitability) . The lines 
indicate the political boundaries of the states and its respective mu- 
nicipalities (GO: state of Goias; MS: state of Mato Grosso do Sul; MT: 
state of Mato Grosso; PR: state of Parana; SP: state of Sao Paulo). The 
records used to calibrate the models are represented by squares, while 
the records used to test the models are represented by triangles. 



A large area of overlap between the models of Lu. 
longipalpis and VL cases was detected. When consider- 
ing all values of climatic suitability of the VL and Lu. 
longipalpis models, 100% of MS territory was suitable 
for the occurrence of vector and disease (Fig. 1). How- 
ever, the northern and northwestern areas of MS were 
suitable for the occurrence of cases and did not show 
high climatic suitability for the occurrence of Lu. lon- 
gipalpis. By superimposing the models of VL cases and 
Lu. longipalpis, considering the higher values of climat- 
ic suitability (> 0.60), areas climatically more favourable 
for the occurrence of Lu. longipalpis and VL cases were 
observed, representing 4\% of MS territory (Fig. 2). This 
overlapping area included most municipalities of intense 
VL transmission (e.g., Aquidauana, Anastacio, Campo 
Grande, Tres Lagoas) and several vulnerable municipali- 
ties for VL transmission in southeast MS (Fig. 2). 

DISCUSSION 

Predictive models indicated that almost the entire 
territory of MS has the potential for VL occurrence. 
In areas more suitable for the occurrence of VL cases 
and less suitable for the occurrence of Lu. longipalpis, 
transmission is possible by Lu. cruzi (Santos et al. 1998). 
This scenario should occur in Corumba (northwest MS) 
and three other municipalities in northern MS (e.g., Sao 
Gabriel do Oeste, Camapua and Cassilandia) where Lu. 
cruzi was detected. The models also indicated that in 
southern areas of MS, there is high climatic suitability 
for the occurrence of Lu. longipalpis, although the spe- 
cies has only been recorded in one municipality based on 




1: Aquidauana (136) 
2: Miranda (27) 
3: Bodoquena (9) 
4: Porto Murlinho (2) 
5: Bonito (28) 
6: Anaat^ (80) 
7: Nloaque (7) 

6: Guia Lopes da Laguna (10) 
9: Jardim (33) 
10:BelaVista(15) 
11: Dols Irmaos do Buriti (14) 
12: SkJroiandia (9) 
13: r^aracaju (11) 
14: Rio Brilfiante (3) 
15: Itapora (0) 
16: Dourados (4) 
17:Caarapd(0) 
18:Juli(0) 
19: Navirai (0) 
20: Amambai (0) 
21 : Nova Alvorada do Sul (0) 
22: Campo Grande (1462) 
23:Rlt>a5doRloPaido(18) 
24: Nova Andradina (0) 
25: Angelica (0) 
26: ivlnhema (0) 
27: Novo Hoilzonte do Sul (0) 
29: Jatel (0) 
29: Taquarusau (0) 
30: Batalpora (0) 
31:Anaunlindia(0) 
32: Bataguassu (0) 
33: Santa Rita do Pardo (0) 
34: Trds Lagoaa (293) 
35: Selvlrta (5) 

Fig 2: overlapping ecological niche models of Lutzomyia longipalpis 
and human cases of visceral leishmaniasis (VL) in Mato Grosso do Sul 
(MS). The overlapping area is shown in red. The white areas indicate 
the absence or predicted climatic suitability smaller values (values < 
0.60). The lines indicate the political boundaries of the states and its 
respective municipalities (GO: state of Goias; MS; MT: state of Mato 
Grosso; PR: state of Parana; SP: state of Sao Paulo). The numbers in- 
dicate the name of the municipalities of MS and the numbers between 
parentheses indicate notified human cases of VL in MS (2003-2012). 
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conducted entomological research. This result illustrates 
a practical application of ecological niche modelling by 
indicating suitable municipalities for vector occurrence 
in cases of entomological research with low capture ef- 
fort or even unrealised sampling. 

Lu. longipalpis was the most abundant species in 
the anthropic environment in municipalities of MS 
such as Campo Grande, where an increase of 60 times 
in the density of this species was registered compared 
to the last sandfly evaluation in 1999/2000 (Oliveira 
et al. 2006a). The factors associated with the increase 
in the occurrence and abundance of Lu. longipalpis in 
MS are only partially understood. Oliveira et al. (2012) 
found a significant positive linear correlation between 
the abundance of Lu. longipalpis and the percentage of 
vegetation cover; moreover, these researchers observed 
the highest abundance of sandflies in houses with chick- 
ens. The attractiveness of chickens to sandflies has often 
been described (Lainson & Rangel 2005). 

Spatial analysis of VL in MS indicated that the ex- 
pansion route of the disease has been from west to east, 
coinciding with the routes of three anthropogenic distur- 
bances of the environment, such as the construction of 
a gas pipeline that began in 1998 and a federal highway 
and railroad that were constructed in the early XX centu- 
ry (Antonialli et al. 2007). The authors hypothesise that 
the intense migration of workers from the VL endemic 
region in the city of Corumba explains the new epide- 
miological pattern of the establishment of this disease 
in MS. Our study indicated two areas of high climatic 
suitability for the occurrence of VL cases. The area near 
Aquidauana, Dois Irmaos do Buriti and Anastacio co- 
incides with high predicted VL incidence rates in MS 
based on geostatistical models (Karagiannis-Voules et 
al. 2013). The incidence of VL in southeast MS is still 
low, but the present study indicates that climatic condi- 
tions are favourable for the occurrence of Lu. longipalpis 
in this area. Thus, municipalities near Nova Alvorada do 
Sul, Angelica and Jatei should be monitored. 

The influence of climatic variables on the occurrence 
of sandflies and VL has been observed in other stud- 
ies (Franke et al. 2002, Peterson & Shaw 2003). In this 
study, the variable that most influenced the models was 
temperature seasonality. Quintana et al. (2013) noted that 
precipitation seasonality and precipitation in the warm- 
est quarter were the variables that most contributed to 
the ENM of Lutzomyia neivai and Lutzomyia migonei in 
Argentina. These variables, as well as temperature sea- 
sonality, are indicative of areas with well-defined dry and 
rainy seasons. In Brazil, the majority of VL cases have 
been recorded in northeastern and Central-West Brazil 
(MS/SVS 2013), which present a well-defined dry sea- 
son. Ecological studies of sandflies in this region have 
indicated that the population density of these insects in- 
creases after the end of the rainy season (Sherlock 1996), 
which has also been observed for Lu. longipalpis in MS 
(Oliveira et al. 2003, Almeida et al. 2010a). 

Other variables, such as deforestation and socioeco- 
nomic factors, which have been linked to the occurrence 
of other sandfly species (Zeilhofer et al. 2008), may also 
influence the occurrence of Lu. longipalpis. Environmen- 



tal disturbances, rural-urban migration, increased popu- 
lation density, a higher amount of solid waste production 
without proper destination and the close proximity of 
residents to domestic animals are some of the factors that 
may facilitate the urbanisation of VL in municipalities 
with more than 100,000 inhabitants (MS 2006). 

A limitation of the niche model for VL cases in 
MS was the use of the coordinates of the municipality 
administrative centre. A more refined analysis could 
be performed with the geographical coordinates of all 
confirmed human cases of VL in MS. Unfortunately, 
this information was not available for the present study. 
Future studies based on the specific location where in- 
fection occurred would increase the number of records 
analysed and allow for a better comparison between 
niche models of vector species and the occurrence of 
human cases. Moreover, studies analysing the effects of 
deforestation and demographic and socioeconomic vari- 
ables for the occurrence of Lu. longipalpis and VL in 
MS may improve the predictive models presented in our 
study and aid in understanding the geography of disease 
transmission in the state. Analysis with higher spatial 
resolution, considering niche characteristics over fu- 
ture climate change scenarios (Moo-Llanes et al. 2013) 
and using different spatial analysis (Sarkar et al. 2010, 
Karagiannis-Voules et al. 2013), may refine the predic- 
tive models presented to better assess the risk factors for 
transmission of VL in MS. 

Species distribution models may contribute to the 
planning of measures designed to control vector-borne 
diseases. Modelling enables the development of geo- 
graphic predictions of pathogens transmitted by insect 
vectors, with important implications for human health. 
Accordingly, the methodological approach presented 
here may serve as a model for studies in other Brazilian 
states. ENM of vectors and human cases allows for a bet- 
ter understanding of the geographic distribution of VL in 
MS, which can be applied to the development of future 
strategies for the surveillance and control of this disease. 
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